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A CONSTRUCTION OF INFINITE BROWNIAN PARTICLE SYSTEMS 


MYKHAYLO SHKOLNIKOV 


Abstract. The paper identifies families of quasi-stationary initial conditions for infinite Brownian par¬ 
ticle systems within a large class and provides a construction of the particle systems themselves started 
from such initial conditions. Examples of particle systems falling into our framework include Brownian 
versions of TASEP-like processes such as the diffusive scaling limit of the g-TASEP process. In this con¬ 
text the spacings between consecutive particles form infinite-dimensional versions of the softly reflected 
Brownian motions recently introduced in the finite-dimensional setting by O’Connell and Ortmann and 
are of independent interest. The proof of the main result is based on intertwining relations satisfied by 
the particle systems involved which can be regarded as infinite-dimensional analogues of the suitably 
generalized Burke’s Theorem. 


1. Introduction 

Over the past decade Brownian particle systems have played a crucial role in a variety of problems 
in probability theory. Examples include (among others) the Dyson Brownian motion which is of great 
importance in the study of universal phenomena in random matrix theory (see the survey [7] and the 
references therein), as well as in obtaining the infinite-dimensional dynamics of eigenvalues arising in 
the bulk and at the edge of the spectrum of large random matrices (see |32] for the former and |26] , |27] , 
[28| for the latter); the stochastic version of the dyadic model which helps to understand anomalous 
dissipation of energy in the context of the Euler equation of fluid dynamics (see 12], 0, HI); the systems 
of sticky Brownian particles giving rise to novel stochastic flows of kernels (see m) ; and rank-based and 
volatility-stabilized models of stochastic portfolio theory offering insights into the degree of stability of 
the capital distribution curve in financial markets m, m, [3Qi, m)- Despite these developments a 
unified approach to existence, uniqueness and (quasi-)stationarity in the case of infinitely many particles 
is missing. 

The aim of this paper is to explicitly identify families of quasi-stationary measures for a large class of 
infinite Brownian particle systems, that is, initial conditions for which the joint distribution of spacings 
between consecutive particles does not change over time, and to give a construction of the particle 
systems started from such initial conditions. In addition, a martingale problem characterization is pro¬ 
vided, yielding a powerful tool for establishing convergence to the infinite Brownian particle systems. 
The spacings between consecutive particles in such systems can be viewed as infinite-dimensional ver¬ 
sions of the softly reflected Brownian motions recently introduced in the finite-dimensional setting by 
O’Connell and Ortmann in [23| and are of independent interest. 

Specifically, we consider systems of infinitely many particles on the real line whose positions 
k gN evolve according to the system of stochastic differential equations (SDEs) 

/c+d—1 

(1.1) dXk{t) = dBk{t) + ^ikdt + X U'{Xi{t) - Xi+i{t)) rikdt, keN. 

l=k 

Here {Bk ■ A: G N} is an inhnite system of one-dimensional Brownian motions with quadratic covaria¬ 
tions {Bk,Bi){t) = akit, t > 0 where aki, {k,l) G are real constants; /x^, A: G N is a sequence of real 
numbers with ^ko — k-ko+i = • • • for some A:o G N; U : [0,oo) — )• M is a differentiable function with a 
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locally bounded derivative; and {l,k) G I > k are real constants. Note that the form of (11.11) 
assumes that the interaction between particles is hierarchical, that is, only particles with a higher index 
enter into the dynamics of any given particle, and of finite range, that is, only finitely many particles 
appear in the dynamics of any given particle. Since a shift of a solution to (jl.ip gives another solution 
of (|l.ll) . it may and will be assumed without loss of generality that Xi(0) = 0. 

Our main assumption is on the quantities aki, {k,l) G and rik, {l,k) G N^, I > k and should 
be viewed as the infinite-dimensional version of the skew-symmetry condition in m- For notational 
convenience we use the convention rik = 0 for I > k d and I < k. 

Assumption 1.1. The covariance matrix A = is such that all submatrices {aki)^i^i, K ^ 'N 

are non-degenerate and 


rik — ruf-j,-!) 2 

(1-2) Uki -|- — ak[i+i) = - - - , (A:, f) G N , I > k, 

(1.3) - ai(fc+i) = ^, A:GN. 

We also assume the normalizations a^k + a(fc+i)(fc+i) — a(^k-ei)k ~ = 1, A: G N and rkk = 1, A: G N. 


Remark 1.2. Assumption 1 1.1 1 is satished if, for instance, 2A is the N x N identity matrix and d = 1 (so 
that rik = 0 for all {l,k) G N^, I > k). This simple setting includes Brownian versions of TASEP-like 
particle systems, and we refer to Section [4] below for a detailed discussion of some examples of this type. 


Remark 1.3. With A := (au -h a(k+i)(i+i) - a(k+i)i - and R := {rik - n(fc+i))^i=i it is easy 

to see from Assumption 11.11 that the N x N matrix 2A — is an upper triangular band matrix with 
band width d. In particular, there is a unique upper triangular matrix {2A — i?) whose formal (left 
and right) multiplication with 2A — R gives the N x N identity matrix. Indeed, viewing the matrix 
(2A — R) as a left inverse first, one can fill in its entries from left to right and from top to bottom 
sequentially by solving the associated linear equations. Then, one checks that the resulting matrix is 
also a right inverse of 2A — R. For future reference we set 


(1.4) 


(vi, V2, (2i - R) ^ - fJ.2,R2 - R3,-- •) 


where the superscript “T” denotes transpose. 


In addition, we make mild assumptions on the interaction potential U. 

Assumption 1.4. It holds Zk ;= /j^exp {2U{z) — 2ukz)dz < oo. A: G N and the probability measures 
^ exp (2 U{z) — 2 Vk z') dz, A: G N have finite second moment and Fisher information: 

(1.5) I — exp (2[/(z) — 2z/fc z) (2 f7'(2) — 2z^fc)^ dz < oo, A:GN. 

Jm. 

Remark 1.5. As will become apparent below, the condition Zk < oo. A; G N ensures the existence of 
finite quasi-stationary measures for the process in (II.ip . Moreover, it is easy to see that condition ()1.5p 
is equivalent to 

(1.6) I exp (2 C/(z) — 2 z/fc z) 17'(z)^ dz < oo, A: G N. 

Jr 


The first main result of the paper reads as follows. 

Theorem 1.6. Suppose that Assumptions and \1.4\ are satisfied. Then there exists a unique Mar¬ 
kovian weak solution of the infinite system of SDEs o whose one-dimensional distributions are given 
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by 


(1.7) 



exp (2 U{xk - Xk+i) -2vk {xk - Xk+i)) d(xfc - Xk+i) 


1 

(27roiii)^/^ 



{xi- 

- 


t > 0. 


Picking a constant // € M, subtracting fid from all fik, A; € N, and replacing U in (11.11) by Ufj,{z) := 
U{z) + fiz (clearly, this has no effect on (ll.ljl ) one obtains the following corollary. 

Corollary 1.7. Suppose that Assumvtion \1.1\ holds. Then, for all fi ^ 'R such that satisfies As- 
sumption \1.4[ the measure 

OO ^ 

(1.8) 5o(dxi) — exp (2 U^{xk - Xk+i) - 2 z/fc {xk - Xk+i)) d{xk - Xk+i) 

k=i 

is quasi-stationary for the dynamics of m- 


In the course of the proof of Theorem 11.61 we establish the following result of independent interest. 

Theorem 1.8. Suppose that Assumptions [7771 and \1.4\ are satisfied and consider a weak solution of 
dni) as in Theorem [13 Then, for any fixed K G N, the process {Xk : 1 < k < K) constitutes the 
unique weak solution of the system of SDEs 

K—1 k-\-d—l 

(1.9) dXk{t) =dWk{t) + fikdt+'^U'{Xi{t) - Xi+i{t))rikdt+ ^ n^ikdt, l<k<K 

l=k l=K 

in its own filtration. Here {Wk : 1 < A: < K) is a Brownian motion with covariance matrix {aki)ki=i- 

The statement of Theorem 11.81 can be interpreted as an averaging principle. Indeed, in the natural 
filtration of the Brownian particle system {Xk : A; E N) the process {Xk : 1 < A; < K) simply satisfies 
the first K SDEs in (II.ip . whereas in the smaller filtration generated by {Xk : 1 < A: < K) (that is, 
after averaging out the interactions with the processes Xi, I > K) it solves (11.91) . That is, as a result of 
the averaging the interaction terms with the processes Xi, I > K get replaced by suitable constants in 
the dynamics of the processes Xk, 1 < k < K . 

Theorem ll.8l provides a powerful tool for proving weak convergence to the solution of (II.ip of Theorem 
11.61 Indeed, in order to prove weak convergence of a sequence of processes (H^ : A; E N), M E N to 
the latter, it suffices to show weak convergence of each of the sequences (E^ : 1 < A: < K), M E N 
to the solution of (jl.9p with the same value of K. This, in turn, can be accomplished using standard 
martingale problem techniques (see e.g. [SI Section 4.8]). 

The proofs of Theorems II.61 and II.81 relv heavily on intertwinings of diffusion processes as discussed in 
great generality in [31]. In fact, as the proofs show, for each K gN, the diffusion of (II.9|) is intertwined 
with the weak solution of (HTP described in Theorem fLHl Moreover, any two diffusions of (II.9p with 
different values of K are intertwined with each other as well. 

The remainder of the paper is structured as follows. Section [2] is devoted to the proof of the existence 
part in Theorem 11.61 The proof is based on the consistency of solutions to (II.9p with different values of 
K which is the subject of Proposition 12.11 below. Section [3] deals with the uniqueness part in Theorem 
11.61 The proof of the latter is based on Theorem 11.81 the proof of which, in turn, makes use of the 
entropy approach to time reversal developed by Follmer and Wakolbinger in the series of papers |9], 
[IDj . [llj and is given in Section [3] as well. Finally, Section [4] provides some examples from the area of 
integrable probability falling into the framework of (II.ip . 
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2. Existence 

A key ingredient in the proof of Theorem 11.61 is the following proposition. 

Proposition 2.1. Suvvose that Assumvtions \l.l\ and \1.4\ are satisfied. Then: 

(a) For each K gN, there exists a unique weak solution of (USD started according to 


( 2 . 1 ) 


(5o(dxi) 


K-l 


n 


1 

'2'fc 


exp (2 U{xk - Xk+i) - 2 Vfc {xk - Xk+i)) d{xk - Xk+i) 


and, for t > 0, its one-dimensional distributions are given by 

[2„antyr- n i .xp{2C/fe)-2n.a)dfe.-x^+i), 


(b) The solution in part (a) has the property that, for any I < J < K, the process : 1 < k < J) 

is a weak solution of (II.9p with K = J. 

Proposition 12.11 yields the existence part of Theorem 11.61 Indeed, it shows that the sequence of 
solutions in G N of (II.9j) with initial conditions of the form (j2.ip is consistent. Hence, one can 

define a process X = {X^, fe G N) as the limit of X^^\ A' G N in the sense of the Kolmogorov Extension 
Theorem (see e.g. m Theorem 5.16]). Due to the fact that the ci-algebras in the natural filtration of 
X are generated by cylinder events, it follows that the processes 

pf k~\~d — 1 

Xk{t)-Xk{0)-fikt- U'{Xi{s)-Xi+,{s))rikds, keN 

^0 l=k 

are continuous martingales in that filtration. Moreover, the quadratic covariations of the latter processes 
are given by Oki t, {k, 1) G and, thus, Ai is a weak solution of (jl.9l) by Levy’s characterization of 
Brownian motion (see e.g. m Theorem 3.3.16]). Lastly, X has the Markov property, since, for any 
t >0, the cr-algebras generated by X on [0,t), at t, and on {t,oo) can be generated by cylinder events, 
and each X^^'l has the Markov property (as a usual consequence of weak uniqueness for any initial 
condition given by (|2.2p . see e.g. m Theorem 5.4.20]). 


Next, we give the proof of Proposition 12.II 


Proof of Proposition [27X1 Step 1. We fix a iL G N throughout and start with the proof of part (a) 
in the case that U' is bounded. The existence and uniqueness of a weak solution X^^') to (II.9|) in this 
case is a well-known consequence of Girsanov’s Theorem (see e.g. 13 Section 5.3 B]). To complete the 
proof of part (a) in this case it therefore suffices to show that, for any t > 0, the distribution of the 
random vector (^x[^\t), Aj'^^(t) — X^[t), ..., X^^}^{t) — X^{t)^ admits the density 

:= - — -j - exp (2 C/fe)-2 








INFINITE PARTICLE SYSTEMS 


5 


To this end, we set 


K-l 


K-l 


= 9 X] + 0(fc+i)(z+i) - a(fc+i)i - afc(z+i)) [p-k - ^fc+i) 


A;,Z=1 
X-l iC-l 


5 


Syfc dyi 

K—1 k+d—1 


+ Y.Y1 "'■'(fc+i)) ^ ^ink-Y.^1 ’^hfc+1)) 7^’ 


fc=l l=k 


k=l 

k+d 


k=l l=K 


l=K 


dyk 


d 


dyk' 


Ck = + i a„ ^ + Ml ^ ^ U'{y,) -■« + E «l n, ^ 


1 = 1 


l=K 


K-l 


+ ^ («lfc “ “l(A:+l)) 


dxi 


k=l 


dxidyk 


and aim to verify the Kolmogorov forward equation ^ = (Ck)*iP- Here {Ck)* is the formal adjoint of 
Ck, viewed as an operator acting on twice continuously differentiable functions. 


A direct computation gives 


A=( 


dt 


1 


2t 2aiit‘^ 


2aii 




Next, we note that 


ic.n = (.4<-‘))-^ M ‘ <111 0 _ Ml w .,1 ^ m,m,i |i 


1=1 


dxi 


l=K 


K-l 


+ ^ {(^Ik — Oi(fc+i)) 


k=l 


dxi 

dxidyk' 


At this point, |24l Theorem 3.2] in conjunction with Assumption 11.11 yield 


K-l 


(^(A- 1)^* J| g^p ^2 u{y^) -2vk Vk) = 0, 


k=l 


so that {A^^ = 0- Calculating the remaining terms directly we get 

1 ^ (xi - yit- Ya=i ni ^ ~ ^ ~ ^fc=i ^ 


iCK)*^=[ -^ + 


K-l 


+ E 


ni 


2aiit^ 

xi- yit- Y)k=i '^ki T^k t 


aiit 


1=1 

K-l 


aiit 


xi- yit- Y)t=i rki T^k t , 
+ ni - , - ii 


l=K 


aiit 


E («ifc - ai(fc+i)) (2 U'{yk) - 2 Vk) 


xi- yit- Ya=i ni T^i t 


k=l 


ant 


'•P- 


Using Assumption 11.11 this can be simplified further to 


1 ^ {xi-yit-Yfi=iriivitf ( ^ 

-h iyi + / ^ fki k'k 


2t 


= \-l+ 


2oiit2 


k=l 


xi- yit- Y,f=i ni n t 

ant 


ip 


2t 2ant^ 


2an 


Ip. 
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Hence, ^ = {Ck)*'^ as desired. In view of the established weak uniqueness and [THl Corollary 1.3, 
implication a) ^ c)] it follows that, for any t > 0, the distribution of the process — 

X^\ ..., at time t has density •). 

Step 2. We now turn to the proof of part (a) in the case that U' is unbounded. To establish weak 
existence we consider a probability space which supports a ii'-dimensional Brownian motion (3 with 
zero drift vector, covariance matrix := {aki)^i=i and initial condition of (j2.1|) . Furthermore, we 
define the stopping times 

(2.3) Tn ■= inf |s > 0 ; |/3i(s) — /3/+i(s)| > n for some 1 < I < K — l|, n G N, 

fix a t G (0,oo), and introduce the measures P^, n G N given by their densities with respect to the 
underlying probability measure: 


dPn 


= exp 


rmin(t,r„) ^ 


K-l 


k+d—1 


'0 


l=K 


f 


- A+i(s)) rik + nrikj J^d/3m(s) 

ds ), n G N. 


/c,m=l l=k 

^min(£,T7T,) . k +d— 1 

, fj-k + u'{/3i{s) - f3i+i{s)) rik + ^ I'lrik 

l=k l=K 


l<k<K (A(^))-i 


Here ^ is the inverse of A^^\ and || • ||^^(x))-i stands for the Euclidean norm on associated 

with . Since U' is locally bounded, P„, n G N are well-defined probability measures by the 

Novikov criterion (see e.g. m Proposition 3.5.12]). 


In addition, one computes 


log^l =E® 

L dQ ^ dQ J 

1 pmm{t,Tn) 

= E''”" - / 

L2 ./n 


log 


dPn 


K—1 k+d—1 

fkk +U'{/3i{s) -/3i+i{s)) rik + ^ i^irik 


l=k 


l=K 



2 1 
ds 

VI 

..se 

VI 

(AW)-i J 


, n G N 


where the second equality is a consequence of Girsanov’s Theorem. At this point, for every fixed n G N, 
one can find a differentiable function Un which coincides with U on [—n, n] such that is bounded, 
Xn k ■= /r exp (2 Un{z) — 2 Vk z) dz < oo, 1 < k < K — 1, and 
(2.4) 

I — — exp (2 Un{z) — 2 i>k z^ U^{z)'^ dz < 1 + I — exp (2 U{z) — 2 Uk z:) U'{z)‘^ dz, 1 < k < K — 1. 
Jr 2,n,k Jr Xk 

Due to the weak uniqueness assertion established in Step 1, under P„, the functional 

K—1 k+d—1 

^^k-k + u'{l3i{s) - Pi+i{s)) rik + ^ i^irikY 


f 


l=k 


l=K 


l<k<K (AW)-i 


ds 


has the same law as the functional 
2 j^rmn{t,Tn) K—1 


f 


k+d—1 


l<k<K (A(^))-i 


,kk+Y^ K{x+>’^{s) - X++{s)) rik + ^ rik 

l=k l=K 

where is the weak solution of (|1.9p with U' replaced by U'^ and initial condition 

K-l ^ 

5o(dxi) — exp{2Un{xk-Xk+i)-2uk{xk-Xk+i))d{xk-Xk+i). 

Z^n,k 


ds 


k=l 
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Relying on the formula for the one-dimensional distributions of obtained in Step 1 one can now 

bound the expectation of the latter functional from above by a constant independent of n using 


E 


f‘K{x. 

Jo 


{K),n, 




= t —— exp (2 Un{z) — Ivi z) dz, 1 < I < K — 1 

Jr ^n,l 


and ()2.4I) . 


All in all, it follows that the quantities log , re G N can be bounded above uniformly 

in re, so that re G N is a uniformly integrable martingale under Q (the martingale property can 
be deduced directly from Ito’s formula). Thus, Doob’s Martingale Convergence Theorem yields the 
existence of the limit ^ := lim„_;.oo in the sense, and we can let P be the probability measure 
associated with the density Moreover, the convergence implies that, for every twice continuously 
differentiable function / : —>• M which is bounded together with all its first and second order partial 

derivatives, re G N, si,S 2 G [0,t], and bounded random variable A measurable with respect to the 
u-algebra generated by /3(min(s, r^)), s G [0, si], one has 



f{/3{mm{s2,Tn))) -/(/3(min(si,rn))) - / 

J IT 


min(s2Tn) 

min(si,r„) 


{CkJ) (/3 (min(s, Tn))) ds ) A 


= 0 


where Ck is the generator corresponding to the SDE (II.9p . In other words, under P, the process /3 solves 
the local martingale problem associated with (II.9p on [0, t] and is, hence, a weak solution of (II.9p on 
[0,t] (see e.g. [171 Corollary 5.4.8]). The weak existence on the whole time interval [0,oo) follows from 
the Kolmogorov Extension Theorem and the weak uniqueness on every time interval [0, t] established 
next. 


In order to show weak uniqueness on each time interval [0,t], we pick a weak solution of (II.9p 

on [0,t], redefine the stopping times Tn, re G N by replacing A(s) — /3i+i{s) with X^^^s) — X^^^ls) in 
(12.31) and only considering s G [0, t] there, and note that 


(2.5) 




= lim E^ 
n—>-oo 


K-1 

E' 

l = k 

rrn 


k-\-d—l 


l=K 


) 

2 

ds 

/ l<fc<A' 

(AW)-i 


Tn K — 1 k-\-d—l 

('5)) nk + '^ nk) 

l=k l=K 



2 

ds 

VI 

VI 



by the Monotone Convergence Theorem. Moreover, the latter expectations can be bounded above 
uniformly in re as before, by replacing U with Un and with changing the interval of 

integration from [0,r„] to [0,t], and using the formula for the one-dimensional distributions of 
obtained in Step 1. It follows that the expectation on the left-hand side of (|2.5p is finite. 

At this point, m Theorem 7.5] shows that the law of A^^^ on [0,t] is absolutely continuous with 
respect to the law of the previously introduced Brownian motion /3 on [0, f]. Finally, the formula for 
the associated change of measure in |19[ Theorem 7.6, equation (7.29)] demonstrates that the law of 
aW is uniquely determined. 


Step 3. To prove part (b) one only needs to consider the case J = K — 1, since for all other values of 
1 < J < A one can apply the statement for J = A — 1 repeatedly. Set 


yiK) ^ x(K) _ y 

FAA-1) ^ ^(K) 


(K) 
2 ) 


^2 ) • • ■ 5 ^K-2 ^K-l) 


^{K) _ ^(K) 
{K) 




AK) 


and let A(f), t > 0 and ^(t), t > 0 be the filtrations generated by and yAA 1 )^ respectively. 
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It suffices to show that the process 

(2.6) M^(t) = - / {CK-if){Y^^’^-^\s)) ds, t>0 

Jo 

is a martingale with respect to G{t), t > 0 for any twice continuously differentiable function / : —)• 

M which is bounded together with all its first and second order partial derivatives. To this end, we 
employ the martingale property of the process 

_ ;(y(i^,i^-i)(o)) _ f\cKf){Y(^){s)) ds, t > 0 

Jo 

with respect to t > 0 (due to Ito’s formula) to compute 


= f{Y^^’^-^Hs)) +e\e\ [\cKf){Y^^\u))du 


ns) 


g{s) 


= f{Y^^'^-^Hs))+E / {CK-if){Y^^’^-^\u)) du 


g{s) 


+E 


{{CK-CK-i)f){Y^'"\u))^u 


g{s) 


for any t > s > 0. Therefore we need to show that the last summand in the latter expression vanishes. 
Moreover, due to Fubini’s Theorem, it is sufficient to prove 


(2.7) 


E 


{iCK-CK-i)f){Y^^\u)) g{s) 


= 0 , 


u> s>d. 


We explain first how (|2.7p can be obtained from the following claim the proof of which we defer to 
Step 4 below. 

Claim 1. For any u > 0, the cr-algebra g{u) and the random variable Y^\u) are independent. 

Given the claim we find 

E\{CK-CK-i)f){Y^'^\u)) g{s) 


g{s) 

= E 

E 

\{CK-CK-i)f){Y^'^\n)) 

g{u) 

g{s) 




= E 

E 

\{CK-CK-i)f)iY^'^\n)) 


)_ 

G{s) 


Hence, to obtain (|2.7p it suffices to prove 


(2.8) E 

At this point, we use formula (j2.2p to write 
E 




= 0, u>d. 


X ^ ^ inK-i)i - nK-m+i)) 


1=1 


■ {U'{z) — vk) exp(2 U{z) — 2 vk z) dz. 

Since the latter integrand is a derivative in z of a function vanishing at positive and negative infinity, 
its integral equals to zero. All in all, it follows that the process of (j2.6h is a martingale with respect 
to ^(t), t > 0 for any function / as described in the paragraph after (12.6p . 


Step 4. To finish the proof it remains to establish Claim 1. To this end, it suffices to consider the time- 
reversed process Y^^’^~^\s) := — s), s G [0, w] in the filtration F{s), s G [0,u] generated 

by Y^^\u — s), s G [0,u] and to verify that its dynamics does not depend on Y^\u). 







































INFINITE PARTICLE SYSTEMS 


9 


Recall from Steps 1 and 2 that there exists an equivalent probability measure Q on the underlying 
probability space under which s € [0,rt] is a Brownian motion with zero drift vector, 

covariance matrix with entries 


hi ■■= 


and initial condition 
(2.9) 


aki 

afc(z-i) - CLki 
<^{k-i)i - o-ki 
®(fc-l)(Z-l) + 0,kl 

K-2 


«fc(Z-l) — <^[k-l)l 


if k = 1 = 1 

if k = 1, I > 1 

if k > 1, I = 1 

if A; > 1, / > 1 


^o(dxi) ^ exp {2U{yk) -2vk Vk) dyfc 


k=l 


in the filtration s G [0,n]. Consequently, under Q and in the filtration F{s), s G [0,rt] the time- 

reversed process s G [0,n] is a Brownian motion with zero drift vector and covariance 

matrix which is conditioned on having the distribution (12.9p at time u. More specifically, 

F('^’^“^)(s), s G [0,u] solves the SDE 


( 2 . 10 ) 

dy(^’^-i)(s) = dw{s) + 0(-^-i) log 


exp 


{y-z)'^{e^^ jy-z) 

2{u — s) 

5o(d2:i) exp (217(2:;) - 2l/; Zi) dz; 


K-l 


1=2 


ds 

=: d5}(s) + 7 ( 2 )(s,y(^’-^-i)(s))ds 

where in is a (K — l)-dimensional Brownian motion with zero drift vector and covariance matrix 0(^“^) 
under Q and in the filtration J-{s), s G [0,tt] (take e.g. both X and B to be in m Theorem 

!])• 

Note that the drift function 7^^^(s,y) can be rewritten as 


u — s 


+ 0(i^-i) Vj^log 


exp 


yT(^Q{K 1 )^ 


u — s 


exp 


’( 0 (^- 1 )) 


-1 


2{u — s) 

(5o(d2;i) exp ( 21 /( 2 :;) - 2i// 2 ;) da:; 


K-l 


1=2 


and that the second summand in the latter expression is given by a linear transformation of the gradient 
of the logarithmic moment generating function of a hnite measure with sub-Gaussian tails. Hence, 
'y^‘^\s,y) is of uniform sub-linear growth in y on compact subsets of s G [0,n). It now follows from 
m Corollary 3.5.16] that the law of Y^^’^ ^H®)) ^ ^ [0)^) under Q in the filtration J-{s), s G [0, u) 
is locally absolutely continuous with respect to the law of w under Q in the filtration X{s), s G [0, u) 
with the corresponding densities being given by 


( 2 . 11 ) 


exp 


s G [0, u). 


On the other hand, one can employ the key observation made in |10l proof of Lemma 3.1]: since the 
law of the process Y^^’^~^\s), s G [0, rt] under the original probability measure P has hnite relative 
entropy with respect to the law of the Brownian motion with the same initial condition, zero drift vector 
and covariance matrix 0(^“^) (see the estimates on the quantity of (j2.5p above), the same is true for 
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the laws of their time reversals. In other words, the law of the process s € [ 0 ,rt] under P 

must have finite relative entropy with respect to the law of the same process under Q. In particular, 
we conclude that the law of the process s G [ 0 , u) under P in the filtration Y{s), s G [ 0 , u) 

is locally absolutely continuous with respect to the law of w under Q in the filtration T{s), s G [ 0 ,u). 
By [T9l Theorem 7.11] it follows that the process Y^^’^~^\s), s G [0,u) satisfies a SDE of the form 

( 2 . 12 ) = 7 (s) ds + d 6 (s), 0 <s<u 

under P and in the filtration s G [0,u), and [191 Theorems 7.5, 7.6] show that the densities 

associated with the latter absolute continuity relation are given by 

(2.13) [0,«). 

Here 7 and b are a non-anticipative process and a Brownian motion with zero drift vector and covariance 
matrix respectively, in the filtration Y{s), s G [ 0 , u). 


To identify the process 7 (and, thus, to complete the proof of Claim 1) one can use the decomposition 
7(5) =: 7^^^(s) + s G [0,u) and combine (|2.11l) . (|2.13p and Girsanov’s Theorem in 

the form of m Theorem 3.5.1] to conclude that the densities associated with the change of measure 
from (I2.10p to (|2.12p are given by 

^^P(/ sG[0,u). 

Hence, the finiteness of the corresponding relative entropy together with the simple localization argu¬ 
ment in 13 proof of Lemma 2.6] imply 




[/' 


|7^^^(9)||(0{x-i))-idg 


< 00, 


s G [ 0 , u) 


and, thus. 


E^ 


[ ||7^^^(9)||^dg 

. Jo 


<00, s G [ 0 , u) 


where jj • jj is the standard Euclidean norm on This estimate can be improved further to 


(2.14) 


E^ 


/ lb(9)ll 

Jo 


dq 


< 00, 


s G [ 0 , u) 


by using the decomposition 7 (g) = ^\q))j recalling the uniform sub-linear growth 

of 7*^^)(g, y(^’^“^)(( 7 )) in Y^^’^~^\q) on compact subsets of g G [ 0 ,tt), and combining part (a) of the 
proposition with the finite second moment condition in Assumption 11.41 


In view of the estimate (I2.14p one can now apply [101 Proposition 2.5] (more precisely its straight¬ 
forward extension to higher dimensions) to obtain the representation 

7 (s) = lim E^+ e) - y(^’-^-^)(s)|.f(s)] = lim f -f{q) dq 

e^-O €4-0 J s —e 

for Lebesgue almost every s G [0, u) where the limits should be understood in the sense. This 
representation can be simplihed further to 

(2.15) 7 ( 5 ) = lim -|-e) — — s)] = lim f 7 (g) dg 

€4-0 €4-0 j g —£ 

by noting that (and, hence, also Y^^\u — s), s G [0, u)) is a Markov process. Indeed, the proof of 
weak uniqueness in part (a) of the proposition can be repeated word for word for an initial condition as 
in (I22D, and this yields the Markov property of by a standard argument (see e.g. m Theorem 
5.4.20]). 
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Now, for any twice continuously differentiable function / ; 
on the one hand, 


with compact support one has, 


(2.16) lim + e) - f {Y^^\u - s)) = E^[ 7 (s) f(Y^^\u - s))] 

for Lebesgue almost every s £ [0,u) by virtue of (I2.15p . On the other hand, writing for the 

drift coefficients of under P and in the filtration t > 0, one computes 


lim e"i E’^ 
etO 


(y(k,k-i)^^ + e) - - s)) 


= - lim e"^ E’^ 
40 

+ lim E’P 
etO 


y(i^,Lf-i)(^ - s) f{Y^^\u - s)) - - s - e) f{Y^^\u - s - e)) 




- s)) - - s - 


= -E^ 


-E*" 


+ E'^ 


y(i^,K-i)(^ - s) {CkI) {Y^^\u - s))j - E^ - s) - s) 




|^{y<«)(u-4ei7 + X:£(rW(.-,))0 


Y(K,K- 1 )^^ _ {y(^)(u - s)) 


{l+l)k J 


l<k<K-l. 


for Lebesgue almost every s £ [0, u) using Ito’s formula and forward in time versions of (j2.15p . By can¬ 
celling out the two identical expectations and integrating by parts one can simplify the latter expression 
to 


(2.17) 

e' 


V V Oll(u-s) 


0 


(K) 

Ik 


K 


1=2 


+ Y^(2U'{Y^^^\u-s))-2u^@ 


'ik 


l<k<K-l 


Since 7 ( 5 ) is measurable with respect to the u-algebra generated by Y^^\u — s), a comparison of (I2.16P 
and p2.17p allows to identify 7 ( 5 ) for Lebesgue almost every s £ [0, u). Moreover, Assumption 11.11 implies 
that, for such s, the drift coefficient 7 ( 5 ) is a deterministic function of Y^^’^~^\u — s) = Y^^’^~^\s) 
only. Finally, using this and repeating the proof of weak uniqueness in Step 2 mutatis mutandis for the 
SDE resulting from (12.121) . one concludes that the dynamics of the process Y^^’^~^\s), s £ [0, u] in 
the filtration Y{s), s £ [0, u] does not depend on Y^\u). Claim 1 readily follows. □ 


3. Uniqueness 

We first explain how Theorem 1 1.81 gives the uniqueness part of Theorem 11.61 and, thus, completes its 
proof. To this end, we consider two weak solutions X and X of pi.II) as described in Theorem 11.61 By 
Theorem Mi for any K £ N, the processes {Xj. ■ 1 < k < K) and {Xj. : 1 < A: < AT) are both weak 
solutions to pi.9p . Therefore, by Proposition [2T] (a) their laws must be the same. It follows that X and 
X have the same law. 

Next, we prove Theorem ll.81 The proof is similar to the proof of Proposition l2.ll but requires certain 
modifications and is therefore presented in full detail for the convenience of the reader. 

Proof of Theorem 11.81 Step 1. We start by fixing a AT £ N, putting 


:= (Ai,Xi-A 2 ,X 2 -A 3 , ...) and := (Ai, Xi - A 2 , ..., 
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and letting t > 0 and t > 0 be the filtrations generated by and 

respectively. Moreover, we formally define the operator 


C 


oo 


1 ^ q 2 d ^ d 

^ ^ + E v'bi)rn 

^ k=l 1=1 


^ cx:> 

+ 9 {cLki + a(fe+i)(;+i) — a(fc+i)i — afe(i+i)) 




kl=l 


dyk dyi 


dxi 

+ ^ {y-k — yk+i) 


k=l 


d 

dyk 


OO OO 

k=ll=k 


'^lik+l)) 


d 

dyk 


which should be understood as acting on functions / : —)■ M depending only on finitely many 

coordinates of (xi, yi, y^, • • •) G and being twice continuously differentiable with respect to these 
coordinates. 


Recalling the notation Ck introduced in Step 1 of the proof of Proposition 12.II we note that, in order 
to establish the theorem, it suffices to show that the process 

(3.1) - /(y(“’^)(o)) - f {CKf){Y^°^'^\s)) ds, t > 0 

Jo 

is a martingale with respect to t > 0 for any twice continuously differentiable function 

/ : —>■ M which is bounded together with all its first and second order derivatives. To this end, we 

use the martingale property of the process 

/(y(“’^)(t)) -/(y(-’^)(o)) - A£oo/)(y(°°)(s))ds, t>o 

Jo 

with respect to t >0 (due to Ito’s formula) to hnd 


E [/ (y(°°’^) (t)) I (s)] = E [E [/ (y (t)) (s)] | (s)] 


= f{Y^°^’^\s))+E f {CKf){Y^°°’^\u))duJ^^°°’^\s 

_ J S 


jr(oo,i^)(^) 


+E 


{{C^-CK)f){Y^°^\u))du 


jr(oo,K)(^) 


for any t > s > 0. It therefore remains to show that the last summand in the latter expression is equal 
to zero. In addition, thanks to Fubini’s Theorem this statement can be further reduced to 


(3.2) 


E 


{{Coo-CK)f){Y^°^\u)) y(°°’'")(s) 


= 0, tt > s > 0. 


We first establish (|3.2I) assuming the following claim and then supply the proof of the claim in Step 
2 below. 

Claim 2. For any u > 0, the cr-algebra and the random vector (Xj^(u) — 

X/c+i(u), Xjc+i(u) — Y/c+ 2 (u), ■ ■ ■) are independent. 


Given the claim we compute 


'((£oo-YK)f)(Y(‘^H^)) 

^(oo,i^)(g) 

= E 

E 

'((£oo-YK)f)(Y(‘^H^)) 


jr(oo,i^)(g) 



= E 

E 

\(£oo-YK)f)(Y(^Hu)) 

y(oo,i^)(^) 

jr(oo,i^)(^) 


































INFINITE PARTICLE SYSTEMS 


13 


Hence, to obtain (|3.2p it suffices to prove 


(3.3) E 

At this point, we use formula (11.71) to write 
E 




= 0, ix > 0. 


rn + X]" ^ [nk - n^k+i)) 


K+d-1 

E 

l=K 


k=l 

K+d-1 


iu'{zi)-n) n 


dvk 

1 


exp(2 U (zm) 2 zxjyj ZjYi) d^^r? 


m=K 


Note that each of the d integrands is a first order partial derivative of a function vanishing at infinity, 
so that all d integrals are equal to zero. The theorem now readily follows. 

Step 2. It remains to supply the proof of Claim 2. To this end, we will analyze the time-reversed process 
y(°°’^)(s) := — s), s e [ 0 ,u] in the filtration s G [ 0 , ix] generated by Y^°°\u — s), 

s G [ 0 , ?x] and verify that its dynamics does not depend on (^Xk{u)—Xk+i{u), Xk+i{u)—Xk+ 2 {u), ...). 

Recalling the notations for the covariance matrix {aki)ki=i and || • ||^^(k))-i for the Euclidean 
norm associated with its inverse we find 


(3.4) 


E 


f 


k-\-d—l 

fik+ U'{Xi{s)-Xi+i{s))rik] 

l=k 



2 

ds 

VI 

VI 

(AW)-i 


< oo 


thanks to Fubini’s Theorem, and dLSD. Therefore, according to m Theorem 7.2], the law of 

the process (Ai(s), A 2 ('S), ..., Xk{s)), s G [0, xx] is absolutely continuous with respect to the law of 
a A'-dimensional Brownian motion with zero drift vector, covariance matrix and initial condition 
(Ai( 0), Ai 2 ( 0 ), ..., Ax(0)). Moreover, by [HI Theorem 7.6, equation (7.30)] the relative entropy 
corresponding to this change of measure is given by one half of the quantity in (13.411 and is, hence, 
finite. In other words, the relative entropy of the law of the process Y^°°’^\s), s G [0,xx] with respect 
to the law of a Brownian motion with zero drift vector, covariance matrix 0(^1 (introduced in Step 4 of 
the proof of Proposition 12 . 11 ) and initial condition jg finite. At this point, as in Step 4 of the 

proof of Proposition 12.II one concludes that the law of the time-reversed process Y^°°’^\s), s G [0,xx] 
has finite relative entropy with respect to the law of a Brownian motion on the time interval [0, xx] with 
zero drift vector and covariance matrix conditioned to have the distribution of y(°°A)( 0 ) at time 

u (that is, the process of (|2.I0p with K replaced by {K -|- I)). 

Recall that in Step 4 of the proof of Proposition 12. II it has been established that the law of the process 
of (I2.IUI) is locally absolutely continuous with respect to the law of a Brownian motion b on the time 
interval [0,xi) with zero drift vector and covariance matrix 0^^). Consequently the law of the process 
y(°°’^)(s), s G [0, xx) is also locally absolutely continuous with respect to the law of such a Brownian 
motion, so that y(oo,A:) j^^g^ satisfy an SDE of the form 

(3.5) dY^°°'^\s) =''y{s) ds + db{s), sG[ 0, xx) 

by [191 Theorem 7.II]. Here 7 is a non-anticipative process in the filtration s G [0, xx). More¬ 

over, according to m Theorems 7.5, 7.6], the densities associated with the latter absolute continuity 
relation are given by 


(3.6) 


exp 


’7((?)^(0W) My(-’^)(g)-^^'||7(g)||^eW)-idg)> sG[0,u). 


To determine 7 we consider the decomposition 7 ( 5 ) =: 7 ^^^(s) -|- {s,Y^°°'^\s)), s G [0, xx) where 

7 ^^) is the same as in (12.10]) . but with {K — I) replaced by K. Combining (|3.6j) with the expression for 
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the change of measure densities in (12.lip (with [K — 1) replaced by K) and Girsanov’s Theorem in the 
form of HZl Theorem 3.5.1] we conclude that the densities associated with the change of measure from 
(12.101) (with {K — 1) replaced by K) to (13.51) are given by 


exp 


^db{q) + - || 7 ^^H 9 )||^e(ir))-idgj, s£[0,u). 


Now, the finiteness of the corresponding relative entropy together with the simple localization argument 
in m proof of Lemma 2.6] give 


E 


|7^^^(g)||JeW)-i dg 


< oo, s £ [ 0 , u), 


so that 


E 


[ dg 

Jo 


< oo, s £ [ 0 , u) 


where || • || is the standard Euclidean norm on This estimate yields further 


(3.7) 


E 


IJQ 


||7(g)|f dg 


< oo, s £ [ 0 , u) 


thanks to 7 (g) = 7 ^^^(g) + 7 *'^^(g,T(°°’^)(g)), the uniform sub-linear growth of in 

F(°°’^)(g) on compact subsets of g G [0,u) (see the paragraph following (12.101) 1. (|1.7I) . and the finite 
second moment condition in Assumption 11.41 

In view of the estimate (13.71) one can now apply cni Proposition 2.5] (more precisely its straightfor¬ 
ward extension to the multidimensional setting at hand) to obtain the representation 

7 (s) = lim E[y(°°’^)(s + e)-y('^’^)(s)|.F(°°)(s)] = lim e-^ f 7 (g)dg 

e4^0 e4'0 Jg —g 

for Lebesgue almost every s £ [0, u) where the limits should be understood in the sense. In view of 
the assumed Markov property of X (and, hence, of yl°°l) the latter representation reduces to 

(3.8) 7 (s) = lim E^ [f (s + e) - T(s) | T {u - s)] = lim f 7 (g) dg. 

c-i-0 J s —£ 


Lastly, pick an L > AT -|- d and a function / : —)• M depending only in the hrst L coordinates 

which is twice continuously differentiable and has compact support with respect to these coordinates. 
Then, on the one hand, 


(3.9) 


lim 

etO 


-1 


E 


(y(oo,K)(^ + e) - - s))j = E[ 7 (s) f{Y^°°\u - s))] 


for Lebesgue almost every s £ [0,tt) by virtue of (|3.8p . On the other hand, with the notations 
for the drift coefficient of with respect to the filtration t > 0 and 


£ 


oo,L 


9 “11 + E(“1^ - “i(fc+i)) TwiV + ^1 TlT: + ^ 


L-1 


92 


d 


dx1 


k=l 


dxidyk 


^x^ 


ni 


1=1 


d 

dxi 


^ L-l q2 

+ 2 X] + «(fc+i)(z+i) - «(fc+i)« - ak{i+i)) Q Q Afc+i) 

k,l=l k=l 


d 

dVk 


L—1 00 

k=ll=k 


'f'lik+l)) 


d 

dVk ’ 


















INFINITE PARTICLE SYSTEMS 


15 


one computes 


lim e ^ E 
eto 


(y(oo,i^)(^ + e) - - s)) 


= - lim e-i E - s) f {u - s)) - - s - e) f {Y^^^u - s - e)) 


etO 


+ lim e"^ E - s - e) 


°^\u-s)) - f{Y^°^\u-s-^ 


= -E 

- E 
+ E 


y(oo,i^)(^ - s) (£ 00 ,l/) (y^°°)(M - s))l - E[/(y(°°)(M - s) - s) 


L-1 


- *)) eSt' + y ^ - *)) 


l<k<K. 




for Lebesgue almost every s G [0, u). Here we have used Ito’s formula and forward in time versions of 
(|3.8j) . By cancelling out the two identical expectations and integrating by parts one can simplify the 
latter expression to 


(3.10) 
E 


/(y(-)(. -S))(- -s) + (- + 

V V aii(u-s) 


0 


(L) 

Ik 


1=2 


+ ^ (21/' (y/“) (u - s)) - 2 ) 0 


hk 


l<k<K 


Since 7(5) is measurable with respect to the u-algebra generated by Y^°°\u — s), a comparison of 
(j3.9l) and (13.101) allows to determine 7(5) for Lebesgue almost every s G [0,rt). Using Assumption 
11.11 one can deduce further that, for such s, the drift coefficient 7(5) is a deterministic function of 
y(oo,iC)(^ - s) = Y^°°’^\s) only. The proof of weak uniqueness in part (a) of Proposition 12.11 can be 
now repeated mutatis mutandis to show weak uniqueness for the SDE resulting from (13.51) . and one 
concludes that the dynamics of the process Y^°°’^\s), s G [0, u] in the filtration s G [0, n] does 

not depend on (^Xk{u) — Xk+i{u), Xk+i{u) — Xx+ 2 (u), ...). This gives Claim 2. □ 


4. Examples 


To illustrate the wide applicability of our results we give two examples of infinite Brownian particle 
systems arising in the area of integrable probability to which they apply. 


Example 4.1 (/3 analogues of Brownian TASEP). Fix a /3 > 0 and consider the /3 corners process which 
is the probability measure on the Gelfand-Tseitlin cone 

{x = (x{: l<i<j<N)e m^(^+i)/2 : < x^l, l<i<j<N-l] 

with density proportional to 


N 




(x. 


N\2 


1<2<J'<A^ 


2 = 1 


N 

n n ' 

j = l 


xl — xi 




j i+1 

nni 

L=l l' = 1 


X1 — X' 


i+1 


|/ 3/2 


-1 


The interest in this measure comes from the fact that, for /3 = 1, 2, and 4, it describes the joint 
distribution of the eigenvalues of the lxl,2x2,...,yxy top left corners of a random matrix from 
the Gaussian orthogonal ensemble (GOE), Gaussian unitary ensemble (GUE), and Gaussian symplectic 
ensemble (GSE), respectively (see e.g. [22] and [T2| for more details). In the papers [33], [12] dynamic 
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versions of the /3 corners process were introduced for /3 = 2 and /? > 2 , respectively, with the latter 
being given by the diffusion process 


dRlit) =dwl{t)-Y^ 


(3/2-1 


■dt + ^ 


(3/2-1 




dt, 1 < i < j < N 


on the Gelftand-Tseitlin cone. Here wj, 1 < i < j < N are independent standard Brownian motions. 

The study of the extremal particles R^, R^Zi, ... is of particular interest due to their relations with 
objects in the Kardar-Parisi-Zhang (KPZ) universality class. Indeed, for (3 = 2, the extremal particles 
form a Brownian analogue of the totally asymmetric simple exclusion process (see [33] for more details) 
which is known to belong to the KPZ universality class (see e.g. [ 6 ] and the references therein). For 
(3 > 2 the situation is more intricate, but one can show that, at least for /3 > 4, the process of extremal 
particles R^, ■ ■ - j seen from the point of view of the particle R^, converges (after the appropriate 

rescaling) in the limit ^ oo to the spacings in the infinite Brownian particle system 

(4.1) = + ..N 

endowed with the quasi-stationary initial condition 

oo ^ 1^ j(2^ 1^ 

(4.2) <5o(dxi) t{(3/2) ( 2 ) ~ exp | {xk - x^+i)) d{xk - Xk+i). 

Hereby, r(-) is the Gamma function. We refer to |13j for more details and note that in (BID time is 
slowed down by a factor of two compared to the setting there to concur with the normalizations of 
Assumption 11.11 


One can put the system ()4.1jl into the framework of dHD by setting /xfc = /r. A: G N for an arbitrary 
/r > 0, U{z) = (4 ~ 5 ) log^ — iJiz, and rik = lk=h {k,l) £ From Remark 11.21 it follows directly 
that Assumption 1 1.1 1 is satisfied. Moreover, the density l(o,oo)(' 2 ^) can be normalized to the 

probability density of the appropriate Gamma distribution (with /i = f corresponding to ()4.2p l. the 
latter has a finite second moment, and its Fisher information is finite iff 

r ^dz<oo, 

Jo 

that is, iff /3 > 4. It follows that, for (3 > 4, the results of Theorems 11.61 and 11.81 apply to the system 
(j4.ip and, in particular, identify all measures 

(4.3) 5o(dxi) r{(3/2) ~ - a^fc+i)) d(xfc - Xk+i), /r > 0 

as quasi-stationary for that system. 


Example 4.2 (O’Connell-Yor semi-discrete polymer Sz Brownian g-TASEP). Fix an G N and a t > 0, 
and imagine a right and up polymer path in the plane connecting (0,0) to (t, N) by following the line 
M X {0} up to the point (ti,0) for some 0 < < t, then moving up to (ti, 1) and following the line 

M X {1} up to the point {t 2 , 1) for some ti <t 2 < t, then moving up to {t 2 , 2) etc. until the point {t, N) 
is reached. In the O’Connell-Yor semi-discrete polymer model (see [25|, m, m) every such path is 
assigned a weight proportional to 

exp -h {bi{t 2 ) - bi{ti)) (6Ar(t) - bN{tN))) 

where bo, bi, ..., b]\f are independent standard Brownian motions. The logarithmic partition function 
associated with this measure is then given by 

ZN{t) := log f exp (bo(ti) + (bi(t 2 ) - bi(ti)) + ... + (bN{t) - bN{tN))) dti dt 2 ... dtN. 

J0<t\<t2<-<tM<t 
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As one lets N £ N and t > 0 vary, one finds (see m) that each process (Zo(t), , Z]y{t)^, 

t > 0 is a diffusion satisfying 

dZfc(t) = d 6 fc(t) + e^'=-iW-^'“Wdt, fe = l, 2, iV, 

dZo(t) = d6o(t). 

The same diffusion also appears in the context of the Whittaker growth model (see [5l Definition 4.1.26]) 
and provides the scaling limit of the g-TASEP process introducted in [5l Section 3.3.2] (see [5l Remark 
4.1.28] for more details). It can be also viewed as a series of queues in tandem as explained in [25]. A 
natural infinite-dimensional analogue of (|4.4I) is the system of SDEs 

(4.5) dXkit) = ^ dbk(t) + ^ dt, A; = 1, 2, ... 

which, in particular, has an interpretation as an infinite series of queues in tandem similar to the one for 
the system (|4.4p . As in the previous example the slowdown by a factor of two is solely for consistency 
with Assumption ll.il 

The system (14.51) falls into the framework of dni) with fi/c = k £ N for an arbitrary > 0, 
[/(z) = —\{^z + e~^), and rik = lk=i, {k,l) £ N^. According to Remark 11.21 the system (14.5p satisfies 
Assumption ll.il In addition, the density exp ( — /j,z — e~^) can be normalized to a probability measure 
which has finite second moment and Fisher information, so that Assumption 11.41 is fulfilled. Thus, 
Theorems 11.61 and 11.81 apply to the system (14.5p and, in particular, show that all measures of the form 

(4.6) 5o(dxi) exp - fi{xk - Xk+i) - d{xk - Xk+i), fJ. > 0 

are quasi-stationary for that system. As before r(-) stands for the Gamma function. 
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